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Abstract 



The /^-differential inclusive production cross section of the prompt charm-strange meson D+ in the 
rapidity range \y\ < 0.5 was measured in proton-proton collisions at ^fs = 7 TeV at the LHC using 
the ALICE detector. The analysis was performed on a data sample of 2.98 x 10 8 events collected with 
a minimum-bias trigger. The corresponding integrated luminosity is L; nt = 4.8 nb _1 . Reconstructing 
the decay D+ — > 07T + , with — > K~K + , and its charge conjugate, about 480 mesons were 
counted, after selection cuts, in the transverse momentum range 2 < pj < 12 GeV/c. The results are 
compared with predictions from models based on perturbative QCD. The ratios of the cross sections 
of four D meson species (namely D°, D + , D* + and D+) were determined both as a function of pi 
and integrated over pj after extrapolating to full pj range, together with the strangeness suppression 
factor in charm fragmentation. The obtained values are found to be compatible within uncertainties 
with those measured by other experiments in e + e~, ep and pp interactions at various centre-of-mass 
energies. 



*See Appendix fAlfor the list of collaboration members 
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1 Introduction 

The measurement of open charm production in proton-proton (pp) collisions at the Large Hadron Col- 
lider (LHC) provides a way to test predictions of quantum chromodynamics (QCD) at the highest avail- 
able collision energies. Charm and beauty production cross sections can be computed in perturbative 
QCD (pQCD) using the factorization approach HUH. In this scheme, cross sections are computed as a 
convolution of three terms: the partem distribution functions of the incoming protons, the partonic hard 
scattering cross section, and the fragmentation process. The partonic hard scattering cross section is 
computed through a perturbative calculation U1I3. while the parton distribution functions and the frag- 
mentation process are parametrized on experimental data. In particular, the fragmentation describes the 
non-perturbative transition of a charm quark to a hadron. It is modeled by a fragmentation function, 
which parametrizes the fraction of quark energy transferred to the produced hadron, and by the fragmen- 
tation fractions, f(c — > D), which describe the probability of a charm quark to hadronize into a particular 
hadron species. 

The production of prompt D°, D + and D* + mesons in pp collisions at y/s = 7 TeV was measured with the 
ALICE detector at two centre-of-mass energies, namely 7 and 2.76 TeV |[3ll4l . Here, 'prompt' indicates D 
mesons produced at the pp interaction point, either directly in the hadronization of the charm quark or in 
strong decays of excited charm resonances. The contribution from weak decays of beauty mesons, which 
give rise to feed-down D mesons displaced from the interaction vertex, was subtracted. The measured pj- 
differential cross sections for prompt D°, D + and D* + are described within uncertainties by theoretical 
predictions based on pQCD at next-to-leading order (e.g. in the general-mass variable-flavour-number 
scheme, GM-VFNS (H) or at fixed order with next-to-leading-log resummation (FONLL (5]). The 
central value of the GM-VFNS predictions for these three mesons lies systematically above the data. On 
the other hand, the data tend to be higher than the central value of the FONLL predictions, as it was 
observed at lower collision energies, namely at the Tevatron (7j[8l, where hadronic decays of D mesons 
were reconstructed, and at RHIC, where measurements of electrons from semileptonic D and B decays 
were performed |[9l[T0l. 

Heavy flavour production in hadronic collisions can be calculated also in the framework of kj -factoriza- 
tion with unintegrated gluon distributions (UGDFs) to account for the transverse momenta of the initial 
partons lfTTT4T4l . Calculations of inclusive production cross section of D mesons based on this approach 
in the leading order (LO) approximation were recently published for LHC energy and compared to 
experimental results 111511161 . 

The measurement of the pj -differential prompt D+ meson production is of particular interest due to its 
strange valence quark content. The D+ production cross section in hadronic collisions was measured at 
lower energies at the Tevatron collider in the transverse momentum (pj) range 8 < pj < 12 GeV/c 0. 
Preliminary results for D+ production at the LHC were reported by the LHCb Collaboration for prompt 
mesons at forward rapidity [171 and by the ATLAS Collaboration at central rapidity |[T8l . The LHCb 
Collaboration also measured the asymmetry between prompt D+ and D~ production in the rapidity 
region 2 < y < 4.5 and for transverse momenta pj > 2 GeV/c, observing a small excess of mesons: 
A P = (ct(D+)-ct(D-))/(ct(D+) + ct(D-)) = (-0.33±0.22±0.10)% ED. Such a particle-antiparticle 
production asymmetry is understood in phenomenological models as due to the effect of the beam 
remnants on the heavy-quark hadronization, see e.g. GUI . 

Charm production has been measured in ep interactions at the HERA collider by the ZEUS ll2ll and 
HI E21 Collaborations, as well as in e + e~ annihilations, at the Z° resonance, by the ALEPH E3l . 
DELPHI EH and OPAL J25l Collaborations, and at centre-of-mass energies of about 10 GeV by the 
CLEO (261 and ARGUS (271 Collaborations. 

As far as theoretical models are concerned, a calculation of the D+ production cross section within 
the FONLL framework is not available, because of the poor knowledge of the parton fragmentation 
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function. The measured data points can be compared with the GM-VFNS prediction that uses meson 
specific fragmentation functions |[28l . 

From the differential production cross section of prompt D°, D + , D* + and D+ mesons, the relative 
production yields of the D meson species can be studied as a function of transverse momentum. A pj 
dependence is expected for these ratios, due to differences in the fragmentation function of the charm 
quark in the four considered meson species, and because of the different contributions from decays of 
higher excited states. In this sense, the measurement of the ratios between the D meson species can 
provide information on the fragmentation functions that can be used in the pQCD models based on the 
factorization approach. The suppression of strange meson production in the charm fragmentation is 
quantified by the strangeness suppression factor, y s , which is computed from the measured D°, D + and 
D+ cross sections extrapolated to full pj range, as defined in Section [6] The values measured at the LHC 
can be compared with those measured for different energies and different colliding systems |[29l . 

Furthermore, the measurement of D+ in pp collisions provides a reference for the studies of charm 
production in heavy-ion collisions. According to QCD calculations on the lattice, under the conditions 
of high energy-density and temperature that are reached in these collisions, the confinement of quarks 
and gluons into hadrons vanishes and a transition to a Quark-Gluon Plasma (QGP) occurs ll30l . Charm 
hadrons are a powerful tool to study the properties of the QCD medium created in these collisions f3l\ - 
l33l . In particular, the meson is sensitive to strangeness production in heavy-ion collisions. Strange 
quarks are abundant in the QGP, resulting in an enhanced production of strange particles with respect 
to pp collisions ll34l437l . Hence, at low momentum, the relative yield of D+ mesons with respect 
to non-strange charm mesons (such as D and D + ) is predicted to be enhanced in nucleus-nucleus 
collisions |T38ti40l . if the dominant mechanism for D meson formation at low/intermediate momenta 
is in-medium hadronization of charm quarks via coalescence with strange quarks ll4TT443Tl . 

In this paper, we report on the measurement of D+ production cross section in pp collisions at y/s = 7 TeV 
with the ALICE detector at the LHC. D+ mesons were reconstructed through their hadronic decay chan- 
nel D+ — > <pn + with a subsequent decay <p — > K~K + . The /^-differential cross section is measured over 
a range of transverse momentum extending from 2 GeV/c up to 12 GeV/c at central rapidity, \y\ < 0.5. 
In Section [2 the detector layout and the data sample are described. This is followed, in Section [3j by 
the description of the D+ meson reconstruction strategy, the selection cuts, and the raw yield extraction 
from the invariant mass distributions. The various corrections applied to obtain the production cross sec- 
tions are illustrated in Section |4] This also includes the estimation of the fraction of promptly produced 
D+ mesons. The various sources of systematic uncertainties are discussed in detail in Section [5] The 
results on the pj -differential cross section compared with pQCD theoretical predictions, the D meson 
production ratios, and the strangeness suppression factor are presented in Section [6] 



2 Detector layout and data collection 

The ALICE detector is described in detail in ll44l . It is composed of a central barrel, a forward muon 
spectrometer, and a set of forward detectors for triggering and event characterization. The detectors of 
the central barrel are located inside a large solenoid magnet that provides a magnetic field B = 0.5 T, 
parallel to the beam line. 

D+ mesons, and their charge conjugates, were reconstructed in the central rapidity region from their 
decays into three charged hadrons (K~K + 7T + ), utilizing the tracking, vertexing and particle identification 
capabilities of the central barrel detectors. 

The trajectories of the decay particles were reconstructed from their hits in the Inner Tracking System 
(ITS) and in the Time Projection Chamber (TPC) detectors in the pseudo-rapidity range |tj| < 0.8. The 
ITS RSI consists of six cylindrical layers of silicon detectors with radii in the range between 3.9 cm 
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and 43.0 cm. The two innermost layers are equipped with Silicon Pixel Detectors (SPD), Silicon Drift 
Detectors (SDD) are used in the two intermediate layers, while the two outermost layers are composed 
of double-sided Silicon Strip Detectors (SSD). The ITS, thanks to the high spatial resolution of the 
reconstructed hits, the low material budget (on average 7.7% of a radiation length for tracks at T] = 0), 
and the small distance of the innermost layer from the beam vacuum tube, provides the capability to 
detect the secondary vertices originating from heavy flavour decays. For this purpose, a key role is 
played by the two layers of SPD detectors, which are located at radial positions of 3.9 and 7.6 cm from 
the beam line and cover the pseudo-rapidity ranges |tj| < 2.0 and |tj| < 1.4, respectively. The TPC ll46l 
provides track reconstruction with up to 159 space points per track in a cylindrical active volume of about 
90 m 3 . The active volume has an inner radius of about 85 cm, an outer radius of about 250 cm, and an 
overall length along the beam direction of 500 cm. 

Particle identification (PID) is provided by the measurement of the specific ionization energy loss, dE/dx, 
in the TPC and of the flight time in the time-of-flight (TOF) detector. The dE/dx samples measured by 
the TPC are reduced, by means of a truncated mean, to a Gaussian distribution with a resolution of 
°d£/dx/ {dE/dx) w 5.5% H6l . The TOF detector is positioned at 370-399 cm from the beam axis and 
covers the full azimuth for the pseudo-rapidity range \r]\ < 0.9. The particle identification is based 
on the difference between the measured time-of-flight and its expected value, computed for each mass 
hypothesis from the track momentum and length. The overall resolution on this difference is about 
160 ps and it includes the detector intrinsic resolution, the contribution from the electronics and the 
calibration, the uncertainty on the start time of the event (i.e. the time of the collision), and the tracking 
and momentum resolution. The start time of the event is defined as the weighted average between the one 
estimated using the particle arrival times at the TOF ATI and the one measured by the TO detector. The 
TO detector is composed of two arrays of Cherenkov counters located on either side of the interaction 
point at +350 cm and —70 cm from the nominal vertex position along the beam-line. In this analysis, 
the time-of-flight measurement provides kaon/pion separation up to a momentum of about 1.5 GeV/c. 

The data sample used for the analysis consists of 298 million minimum-bias (MB) pp collisions at 
y/s = 7 TeV, corresponding to an integrated luminosity L; nt = 4.8 nb _1 , collected during the 2010 
LHC run period. The minimum-bias trigger was based on the information of the SPD and the VZERO 
detectors. The VZERO detector is composed of two arrays of scintillator tiles with full azimuthal 
coverage in the pseudo-rapidity regions 2.8 < r\ < 5.1 and —3.7 < r\ < — 1.7. Minimum-bias collisions 
were triggered by requiring at least one hit in either of the VZERO counters or in the SPD (|t]| < 2), 
in coincidence with the arrival of proton bunches from both directions. This trigger was estimated to 
be sensitive to about 87% of the pp inelastic cross section 1148 11491 . It was verified by means of Monte 
Carlo simulations based on the PYTHIA 6.4.21 event generator i50l (with Perugia-0 tune [51]) that 
the minimum-bias trigger is 100% efficient for events containing D mesons with pj > 1 GeV/c and 
\y\ < 0.5 0. Events were further selected offline to remove the contamination from beam-induced 
background using the timing information from the VZERO and the correlation between the number of 
hits and track segments (tracklets) in the SPD detector. 

During the pp run, the luminosity in the ALICE experiment was limited to 0.6-1.2 x 10 29 cm~ 2 s _1 by 
displacing the beams in the transverse plane by 3.8 times the r.m.s. of their transverse profile, thus keeping 
the probability of collision pile-up below 4% per triggered event. The luminous region, measured from 
the distribution of the reconstructed interaction vertices, had an r.m.s. width of about 4-6 cm along the 
beam direction and 35-50 /im in the transverse plane (the quoted ranges originate from the variations 
of the beam conditions during the data taking). Only events with a vertex found within ±10 cm from 
the centre of the detector along the beam line were used for the analysis. This requirement selects 
a region where the vertex reconstruction efficiency is independent of its position along the beam line 
and it provides almost uniform acceptance for particles within the pseudo-rapidity range |t] | < 0.8 for all 
events in the analyzed sample. Pile-up events were identified by the presence of more than one interaction 
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vertex reconstructed by matching hits in the two SPD layers (tracklets). An event was rejected from the 
analyzed data sample if a second interaction vertex was found, it had at least 3 associated tracklets, and 
it was separated from the first one by more than 8 mm. The remaining undetected pile-up is negligible 
for the analysis described in this paper. 

3 D+ meson reconstruction and selection 

mesons and their antiparticles were reconstructed in the decay chain D+ — > <p7i + (and its charge 
conjugate) followed by — > K~K + . The branching ratio (BR) of the chain D+ — > <p7i + — > K~K + 7T + is 
2.28 ±0.12% ll52l . It should be noted that other D+ meson decay channels can give rise to the same 
KTK+7T+ final state. Among them, those with larger BR are D+ — > K K + and D+ — > f (980)7T + , with 
BR into the KTK+7T+ final state of 2.63±0.13% and 1.16±0.32%, respectively. However, as it will 
be discussed in the following, the selection efficiency for these decay modes is strongly suppressed by 
the cuts applied to select the signal candidates Q, and therefore the measured yield is dominated by the 
D+ — > §k + — > K~K + 7T + decays. 

D+ mesons have a mean proper decay length c% = 150 ± 2 /xm ll52l . which makes it possible to resolve 
their decay vertex from the interaction (primary) vertex. The analysis strategy for the extraction of 
the signal from the large combinatorial background can therefore be based on the reconstruction and 
selection of secondary vertex topologies with significant separation from the primary vertex. 

D+ meson candidates were defined from triplets of tracks with proper charge sign combination. Tracks 
were selected requiring \r[\ < 0.8, pj > 0.4 GeV/c, a minimum of 70 associated space points in the 
TPC, # 2 /ndf < 2 for the track momentum fit in the TPC, and at least 2 associated hits in the ITS, out of 
which at least one has to be in either of the two SPD layers. For tracks that satisfy these TPC and ITS 
selection criteria, the transverse momentum resolution is better than 1% at pj = 1 GeV/c and about 2% 
at pj = 10 GeV/c. The resolution on the track impact parameter (i.e. the distance of closest approach 
of the track to the primary interaction vertex) in the bending plane (r(p) is better than 75 /im for pj > 1 
GeV/c, well reproduced in Monte Carlo simulations 

For each D+ candidate, in order to have an unbiased estimate of the interaction vertex, the event primary 
vertex was recalculated from the reconstructed tracks after excluding the candidate decay tracks. The 
secondary vertex was reconstructed from the decay tracks with the same algorithm used to compute the 
primary vertex [3J. The position resolution on the D+ decay vertices was estimated via Monte Carlo 
simulations to be of the order of 100 /im for each of the three coordinates with little dependence on pj. 
The resolution on the position of the primary vertex depends on the event multiplicity: for the transverse 
coordinates, where the information on the position and spread of the luminous region is used to constrain 
the vertex fit, it ranges from 40 /im in low-multiplicity events to about 10 jltm in events with 40 charged 
particles per unit of rapidity. 

Candidates were then filtered by applying kinematical and topological cuts together with particle iden- 
tification criteria. With the track selection described above, the acceptance in rapidity for D mesons 
drops steeply to zero for |y|>0.5 at low pj and |y|>0.8 at pj>5 GeV/c. A /^-dependent fiducial 
acceptance cut was therefore applied on the D meson rapidity, \y\ < ymipj), where pj is the D+ trans- 
verse momentum. The cut value, yM{pj), increases from 0.5 to 0.8 in the transverse momentum range 
< Pt <5 GeV/c according to a second-order polynomial function and it takes a constant value of 0.8 
for p T > 5 GeV/c. 

The topological selections were tuned to have a large statistical significance of the signal, while keeping 

1 To reduce the combinatorial background, a selection exploiting the mass of the intermediate resonant state was applied. 

t 

Since the width of the <j) peak is narrower than those of the K and the fo(980), the decay channel through the tj> resonance, 
being the one that provides the best discrimination between signal and background, was used in the analysis. 
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Table 1: Measured raw yields (N D > raw ), signal (S) over background (B) and statistical significance ( SA/S + B ) 
for and their antiparticles in the four considered pj intervals. The estimation of the systematic uncertainty on 
the raw yield is described in Section[5] 



pj interval 


N Vf raw ± stat ± syst 


S/B (3a) 


Significance (3a) 


(GeV/c) 








2-4 


125±36±25 


0.12 


3.6 


4-6 


190±35±28 


0.26 


6.3 


6-8 


79±19±12 


0.40 


4.8 


8-12 


85±16±17 


0.58 


5.6 



the selection efficiency as high as possible. It was also checked that background fluctuations were not 
causing a distortion in the signal line shape by verifying that the D+ meson mass and its resolution were 
in agreement with the Particle Data Group (PDG) value (1.969 GeV/c 2 (52 1) and the simulation results, 
respectively. The resulting cut values depend on the transverse momentum of the candidate. 

The candidates were selected according to the decay length and the cosine of the pointing angle, (pointing, 
which is the angle between the reconstructed D meson momentum and the line connecting the primary 
and secondary vertex. The three tracks composing the candidate triplet were required to have small 
distance to the reconstructed decay vertex. In addition, D+ candidates were selected by requiring that 
one of the two pairs of opposite-charged tracks has an invariant mass compatible with the PDG world 
average for the mass (1.019 GeV/c 2 1521 ). To further suppress the combinatorial background, the 
angles Q*{n) and d'(K) were exploited. Q*{n) is the angle between the pion in the KK7T rest frame 
and the KK7T flight line, which is defined by the positions of the primary and secondary vertices in the 
laboratory frame. d'(K) is the angle between one of the kaons and the pion in the KK rest frame. The 
cut values used for the D+ mesons with 2 < pj < 4 GeV/c were: decay length larger than 350 /im, 
cos0 pointing > 0.94, \M™ K _ -M^ DG | < 8 MeV/c 2 , cos0*(7r) < 0.95, and |cos 3 0'(K)| > 0.1. A looser 
selection was applied at higher pj due to the lower combinatorial background, resulting in a selection 
efficiency that increases with increasing pj. 

Particle identification selections, based on the specific energy loss, dE/dx, from the TPC and the time- 
of-flight from the TOF detector, were used to obtain further reduction of the background. Compatibility 
cuts were applied to the difference between the measured signals and those expected for a pion or a 
kaon. A track was considered compatible with the kaon or pion hypothesis if both its dE/dx and time- 
of-flight were within 3a from the expected values, with at least one of them within 2a. Tracks without 
a TOF signal were identified using only the TPC information and requiring a 2a compatibility with 
the expected dE/dx. Candidate triplets were required to have two tracks compatible with the kaon 
hypothesis and one with the pion hypothesis. In addition, since the decay particle with opposite charge 
sign has to be a kaon, a triplet was rejected if the opposite-sign track was not compatible with the kaon 
hypothesis. This particle identification strategy preserves more than 90% of the D+ signal and provides a 
reduction of the combinatorial background under the D+ peak by a factor of 10 in the lowest pj interval 
(2 < ;?t < 4 GeV/c), a factor of 5 in 4 < pj < 6 GeV/c and a factor of 2 at higher transverse momenta. 

For each candidate, two values of invariant mass can be computed, corresponding to the two possible 
assignments of the kaon and pion mass to the two same-sign tracks. Signal candidates with wrong mass 
assignment to the same-sign tracks would give rise to a contribution to the invariant mass distributions 
that could potentially introduce a bias in the measured raw yield of D+ mesons. It was verified, both in 
data and in simulations, that this contribution is reduced to a negligible level by the particle identification 
selection and by the requirement that the invariant mass of the two tracks identified as kaons is compatible 
with the PDG mass. 

The raw signal yields were extracted by fitting the invariant mass distributions in each pj interval as 
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Figure 1: Invariant mass distributions for D+ candidates and charge conjugates in the four considered pj intervals. 
The fit functions described in the text are also shown. The values of mean (jx) and width (a) of the signal peak are 
reported together with the signal counts (S) integrated in ±3c around the centroid of the Gaussian. 



shown in Fig. Q] The fitting function consists of a sum of a Gaussian and an exponential function to 
describe the signal and the background, respectively. For all pj intervals, the invariant mass range 
used for the fit was 1.88 < < 2.16 GeV/c 2 , chosen in order to exclude the region where the 

background shape is affected by D + — > K~K + 7T + decays (BR=0.265% (52)) that give rise to a bump at 
the D + invariant mass (1.870 GeV/c 2 (H). The mean values of the Gaussian functions in all transverse 
momentum intervals were found to be compatible within the uncertainties with the PDG world average 
for the D+ mass. The Gaussian widths are well reproduced in Monte Carlo simulations. The raw 
yield N T>s raw (sum of particles and antiparticles) was defined as the integral of the Gaussian. The 
values of N Ds raw are reported in Table Q] for the different pj intervals, together with the signal-over- 
background (S/B) ratios and the statistical significance, S/\/S + B. For the latter two quantities, signal 
(S) and background (B) were evaluated by integrating the fit functions in ±3a around the centroid of the 
Gaussian. 



4 Corrections 

In order to obtain the pj -differential cross section for prompt (i.e. not coming from weak decays 
of beauty mesons) mesons, the raw yields obtained from the invariant mass analysis (N D > raw ) 
were corrected for the experimental acceptance, the reconstruction and selection efficiency, and for the 
contribution to the D+ measured yield from B meson decay feed-down. The production cross section of 
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Figure 2: Acceptance x efficiency for D+ mesons as a function of pj, for prompt and feed-down D+ mesons (left 
panel) and decays through <p and K intermediate resonant state (right panel). 



prompt D+ mesons was computed as: 



da L 



dpT 



1 1 



prompt 



■N 



raw 



y <yfld 



y|<0.5 



2 Ap T (Acc x e) pr0 mpt ■ BR • L mt 



(1) 



where Apj is the width of the pj interval, Ay (= 2jfld(PT)) is the width of the fiducial rapidity coverage 
(see Section [3]) and BR is the decay branching ratio (2.28% [52 J). The factor / pr0 mpt is the prompt 
fraction of the raw yield; (Acc x £) pr0 mpt is the acceptance times efficiency of promptly produced D+ 
mesons. The efficiency e accounts for vertex reconstruction, track reconstruction and selection, and for 
candidate selection with the topological and particle identification criteria described in Section [3] 
The factor 1/2 accounts for the fact that the measured raw yields are the sum of D+ and DjT, while the 
cross section is given for particles only, neglecting the small particle-antiparticle production asymmetry 
observed by LHCb |fl~9l . The integrated luminosity, L mt = 4.8 nb _1 , was computed from the number 
of analyzed events and the cross section of pp collisions passing the minimum-bias trigger condition 
defined in Section [2 a pp .MB = 62.2 mb 1149 115 31 . The value of ct pp .mb was derived from a van der Meer 
scan ll54l measurement, which has an uncertainty of 3.5%, mainly due to the uncertainties on the beam 
intensities. 

The acceptance and efficiency correction factors were determined using pp collisions simulated with the 
PYTHIA 6.4.21 event generator |[50l with the Perugia-0 tune [51]. Only events containing D mesons 
were transported through the apparatus (using the GEANT3 transport code ll55l ) and reconstructed. The 
luminous region distribution and the conditions (active channels, gain, noise level, and alignment) of all 
the ALICE detectors were included in the simulations, considering also their evolution with time during 
the 2010 LHC run. 



The acceptance-times-efficiency for D+ — > (j)K — > K~K + 7T + decays in the fiducial rapidity range de- 
scribed in Section [3] are shown in the left panel Fig. [2] for prompt and feed-down D+ mesons. The 
acceptance-times-efficiency for the prompt mesons increases from about 1 % in the lowest considered pj 
interval up to 10-15% at high pj. For D+ mesons from B decays, the efficiency is larger by a factor 
1.5-2 (depending on pj) because the decay vertices of the feed-down D mesons are more displaced from 
the primary vertex and, therefore, they are more efficiently selected by the topological cuts. The differ- 
ence between the prompt and feed-down efficiencies decreases with increasing pj, because the applied 
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selections are looser in the higher transverse momentum intervals. The acceptance-times-efficiency for 
prompt D+ mesons obtained without applying the particle identification selection is also shown to sin- 
gle out the PID contribution to the overall efficiency. The used particle identification strategy preserves 
more than 90% of the signal and does not show any significant dependence on D+ meson pj in the range 
considered in this analysis. 

As discussed in Section [3l the decay of the D+ meson into the K~K + 7T + final state occurs via different 

intermediate resonant states. The selection strategy used in this analysis requires that one of the 

opposite-sign pairs of tracks composing the candidate triplet has an invariant mass compatible with the 

meson. The decays D+ — > <p7i + — > K~K + 7T + are therefore preferentially selected by the applied cuts. 

Nevertheless, a fraction of the D+ decaying via another resonant state can pass the selection cuts. In the 

right panel of Fig. [2l the acceptance-times-efficiencies for prompt D+ decaying to KrK + 7T + final state 
— *o 

via a § and a K in the intermediate state are compared. The acceptance-times-efficiency for the decay 
chain D+ — > K K + — > K~K + 7T + is smaller by a factor m 100 with respect to the decay through (j), and it 

is further reduced when applying the PID selection. Indeed, the PID allows the rejection of D+ decaying 
— *o 

via a K that would pass the selection on the invariant mass of the in case of wrong assignment of the 
mass (kaon/pion) to the two same-sign tracks. 

The contribution to the inclusive raw yields due to D+ from B feed-down was subtracted using the 
beauty production cross section from the FONLL calculation HUH, tne ^~ * ®f decay kinematics from 
the EvtGen package Il56l . and the Monte Carlo efficiencies for feed-down D+ mesons. Before running 
the EvtGen decayer, the B admixture cross section predicted by FONLL was split into that of B°, B + , B° 
and Ab by assuming the same pj shape for all hadrons and the production fractions from Il52l . namely 
40.1% of B°, 40.1% of B+, 10.5% of B ( s ' and 9.3% of beauty baryons. The resulting fraction of prompt 
D+ mesons, / prom pt> depends on the pj interval, on the applied selection cuts, and on the parameters used 
in the FONLL calculation for the B meson cross section. It ranges from 0.93 in the lowest transverse 
momentum interval (2 < p T < 4 GeV/c) to 0.87 at high pj (> 6 GeV/c). 

5 Systematic uncertainties 

The systematic uncertainties on the D+ cross section are summarized in Table [2] for the considered pj 
intervals. 

The systematic uncertainty on the yield extraction was defined as the full spread of the D+ yield values 
obtained with different techniques to analyze the invariant mass distributions in each pj interval. The fit 
was repeated in different mass ranges and by varying the function used to describe the background. In 
particular, first and second order polynomials were used instead of an exponential for the background. 
In case of fitting in an extended mass range, a second Gaussian signal was included in the fit function to 
account for the D + — > KrK + 7T + decays. Furthermore, the yield extraction was repeated using a method 
based on bin counting after subtraction of the background estimated from a fit in the mass side bands. 
The resulting uncertainty amounts to 15-20% depending on the pj interval, as detailed in Tabled 

The systematic uncertainty on the tracking efficiency (including the effect of the track selection) was 
evaluated by comparing the probability of track finding in the TPC and track prolongation from the TPC 
to the ITS in the data with those in the simulation, and by varying the track quality selections. The 
estimated uncertainty is 4% per track, which results in 12% for the three-body decay of D+ mesons. 

Another source of systematic uncertainty originates from the residual discrepancies between data and 
simulation for the variables used to select the D+ candidates. The distributions of these variables were 
compared for candidates passing loose topological cuts, i.e. essentially background candidates, and 
found to be well described in the simulation. The effect of the imperfect implementation of the detector 
description in the Monte Carlo simulations was estimated by repeating the analysis with different sets of 
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cuts. The cut values were changed in order to vary the efficiency of signal selection by at least 20% in all 
pj intervals. A systematic uncertainty of 15% was estimated from the spread of the resulting corrected 
yields. Part of this uncertainty is due to residual detector misalignment effects not fully described in the 
simulation. To estimate this contribution, the secondary vertices in the simulation were reconstructed 
also after a track-by-track scaling of the impact parameter residuals with respect to their true value. In 
particular, a scaling factor of 1.08, tuned to reproduce the impact parameter resolution observed in the 
data (see 11310. was used. The resulting variation of the efficiency was found to be 4% in the lowest pj 
interval used in this analysis and less than 1% for pj > 6 GeV/c. This contribution was not included 
explicitly in the systematic uncertainty, because it is already accounted for in the cut variation study. 

Due to the limited statistics, it was not possible to analyze separately D+ and D.T candidates to verify the 
absence of biases coming from a different reconstruction efficiency for tracks with positive and negative 
charge sign not properly described in the simulation H. This check was earned out for other D meson 
species [3j without observing any significant difference between particle and antiparticle. 

The systematic uncertainty induced by a different efficiency for particle identification in data and 
simulation was evaluated by comparing the resulting pj -differential cross section with that obtained 
using a different PID approach based on 3a (instead of 2a) cuts on TPC dE/dx and time-of-flight 
signals, which preserves almost 100% of the signal. In addition, the PID efficiency, was estimated by 
comparing the reduction of signal yield due to the PID selection in data and in simulation, when the same 
topological cuts are applied. Due to the limited statistical significance, this check could be performed in 
data only for D+ candidates integrated over the transverse momentum range 4 < pj < 12 GeV/c. From 
these studies, a systematic uncertainty of 7%, independent of pj, was assigned to the PID selection. 

The contribution to the measured yield from D+ decaying into the KrK + 7T + final state via other resonant 
channels (i.e. not via a meson) was found to be less than 1% due to the much lower selection efficiency, 
as shown in the right panel of Fig. [2] for the case of the decay through a K . The contamination from 
other decay chains (all having smaller branching ratio than the two reported in Fig. [2]) was also found to 
be negligible. 

The effect on the selection efficiency due to the shape of the D+ pj spectrum used in the simulation was 
estimated from the relative difference between the Monte Carlo efficiencies obtained using two different 
pj shapes, namely those from PYTHIA |[50l with Perugia-0 tune BTl and from the FONLL pQCD 
calculation HUH. The resulting contribution to the systematic uncertainty was found to be 3% in the two 
lowest pj intervals, where the selection efficiency is strongly pj dependent, and 2% at higher pj. 

The systematic uncertainty from the subtraction of feed-down D mesons was estimated following the 
same approach as used for D°, D + and D* + mesons 0. The contribution of the FONLL perturbative 
uncertainties was included by varying the heavy-quark masses and the factorization and renormalization 
scales, /ip and /Ir, independently in the ranges 0.5 < \lvjmj < 2, 0.5 < jUr/ott < 2, with the constraint 



0.5 < /If/Mr < 2, where mj = y p T + m\. The mass of the b quark was varied within 4.5 < < 

5 GeV/c 2 . The uncertainty related to the B decay kinematics was estimated from the difference between 
the results obtained using PYTHIA lT50l instead of EvtGen [56] for the particle decays and was found to 
be negligible with respect to the uncertainty on the B meson cross section in FONLL. Furthermore, the 
prompt fraction obtained in each pj interval was compared with the results of a different procedure in 
which the FONLL cross sections for prompt and feed-down D mesons and their respective Monte Carlo 
efficiencies are the input for evaluating the correction factor. Since FONLL does not have a specific 
prediction for D+ mesons, four different approaches were used to compute the pj -differential cross 
section of promptly produced D+. The first two approaches used the FONLL prediction for the generic 
admixture of charm hadrons and that for D* + mesons (the D* + mass being close to that of the D+) 




2 The small particle-antiparticle asymmetry reported by the LHCb Collaboration |19| is negligible in this context. 
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scaled with the fragmentation fractions of charm quarks in the different hadronic species, f(c — > D), 
measured by ALEPH j23l . The other two predictions for prompt D+ were computed using the pj- 
differential cross section of c quarks from FONLL, the fractions f(c — > D) from ALEPH 11231 . and the 
fragmentation functions from l57l . which have one parameter, r. Two definitions were considered for 
the r parameter: i) r = (/md — m c )/mu (mo and m c being the masses of the considered D meson species 
and of the c quark, respectively) as proposed in (571; ii) r = 0.1 for all mesons, as done in FONLL after 
fitting the analytical forms of 11571 to the D* + fragmentation function measured by ALEPH 11581 . The 
D* + mesons produced in the c quark fragmentation were made to decay with PYTHIA and the resulting 
D+ were summed to the primary ones to obtain the prompt yield. For all the four predictions used for 
prompt D+ cross section, the evaluation of / pr0 mpt included the FONLL perturbative uncertainties from 
the variation of the factorization and renormalization scales in the range quoted above and of the c quark 
mass within 1.3 < m c < 1.7 GeV/c 2 . The systematic uncertainty on the B feed-down was defined from 
the envelope of the resulting values of /prompt- The resulting uncertainties in the transverse momentum 
intervals used in this analysis are about , as it can be seen in Table|2] 

Finally, the results have global systematic uncertainties due to the D+ — > (j)7l + — > K~K + 7T + branching 
ratio (5.3% 11521 ) and to the determination of the cross section of pp collisions passing the minimum-bias 
trigger condition (3.5%). 



Table 2: Relative systematic uncertainties for the four considered pj intervals. 





pi interval (GeV / c) 




2-4 


4-6 


6-8 


8-12 


Raw yield extraction 


20% 


15% 


15% 


20% 


Tracking efficiency 


12% 


12% 


12% 


12% 


Topological selection efficiency 


15% 


15% 


15% 


15% 


PID efficiency 


7% 


7% 


7% 


7% 


MC pj shape 


3% 


3% 


2% 


2% 


Other resonant channels 


<1% 


<1% 


<1% 


<1% 


Feed-down from B 


+ 4 % 

-18 10 








Branching ratio 


5.3% 


Normalization 


3.5% 



6 Results 

6.1 /^-differential D+ cross section and D meson ratios 

The inclusive production cross section for prompt D+ mesons in four transverse momentum intervals 
in the range 2 < pj < 12 GeV/c is shown in Fig. [3] As discussed in section [4j the cross section 
reported in Fig. prefers to particles only, being computed as the average of particles and antiparticles 
under the assumption that the production cross section is the same for D+ and D.T . The vertical error 
bars represent the statistical uncertainties, while the systematic uncertainties are shown as boxes around 
the data points. The symbols are positioned horizontally at the centre of each pj interval, with the 
horizontal bars representing the width of the pj interval. In Tabled the numerical values of the prompt 
D+ production cross section are reported together with the average pj of D+ mesons in each transverse 
momentum interval. The (pj) values were obtained from the pj distribution of the candidates in the D+ 
peak region, after subtracting the background contribution estimated from the side bands of the invariant 
mass distribution. The measured differential production cross section is compared to two theoretical 
predictions, namely the GM-VFNS model |6l[59] and the calculations from lfT6ll60l based on the kj- 
factorization approach. 

The GM-VFNS prediction is found to be compatible with the measurements, within the uncertainties. 
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Figure 3: (colour online) ^-differential inclusive cross section for prompt D+ meson production in pp collisions 
at y/s = 7 TeV. The symbols are positioned horizontally at the centre of each pj interval. The horizontal error bars 
represent the pj interval width. The normalization uncertainty (3.5% from the minimum-bias cross section and 
5.3% from the branching ratio uncertainties) is not shown. Theoretical predictions from GM-VFNS (6) and from 
^T-factorization at LO [ 16 , 60] are also shown. 



The central value of the GM-VFNS prediction corresponds to the default values of the renormalization 
(/Ir) and factorization (/ii and /ip for initial- and final-state singularities, respectively) scales, i.e. iir = 

jUi = /ip = mj, where mj = \l p\ + m%, with m c = 1.5 GeV/c 2 . The theoretical uncertainties are 
determined by varying the values of the renormalization and factorization scales by a factor of two 
up and down with the constraint that any ratio of the scale parameters should be smaller than or equal to 
two (H. The central value of the GM-VFNS prediction is higher than the measured point by w 50% in 
the first pj interval, while in the other intervals it agrees with the data within 15%. For D°, D + and 
D* + mesons measured by ALICE at the same pp collision energy Q, the central value of the GM-VFNS 
predictions was found to lie systematically above the data. As mentioned in SectionfjJ predictions for the 
D+ production cross section within the FONLL framework are not available, due to the poor knowledge 
of the fragmentation function for charm-strange mesons. 

The prediction from |[T6ll60l is obtained in the framework of kj -factorization at LO using Kimber- 
Martin-Ryskin (KMR) unintegrated gluon distributions in the proton. The measured D+ cross section is 
described by the upper limit of the theoretical uncertainty band. 

The ratios of the pj -differential cross sections of D + and D* + to that of D°, taken from fj3j, are shown in 
the top panels of Fig. |H In the bottom panels of the same figure, the ratios of the D+ cross section 
to the D° and D + ones are displayed. In the evaluation of the systematic uncertainties on the D 
meson ratios, the sources of correlated and uncorrelated systematic effects were treated separately. In 
particular, the contributions of the yield extraction, cut efficiency and PID selection were considered 
as uncorrelated and summed in quadrature. The systematic uncertainty on the B feed-down subtraction, 
being completely correlated, was estimated from the spread of the cross section ratios obtained by varying 
the factorization and renormalization scales and the heavy quark mass in FONLL coherently for all 
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Table 3: Production cross section in \y\ < 0.5 for prompt mesons in pp collisions at y/s = 7 TeV, in pj intervals. 
The normalization uncertainty (3.5% from the minimum-bias cross section and 5.3% from the branching ratio) is 
not included in the systematic uncertainties reported in the table. The average pj of D+ mesons in each transverse 
momentum interval is also reported. 



pj interval 


(pt) 


da/dpj | y | <0 5 ±stat.±syst 


(GeV/c) 


(GeV/c) 


(Alb GeV^c) 


2-4 


2.7 ±0.4 


19.9±6.1+^ 


4-6 


4.7±0.1 


5.06 ± 1.03+1 j 


6-8 


6.8±0.1 


1.02 ±0.28+^ 


8-12 


9.4±0.1 


0.28 ±0.061°;^ 





P T (GeV/c) p T (GeV/c) 



Figure 4: Ratios of D meson production cross sections as a function of pj. Predictions from FONLL, GM-VFNS 
and PYTHIA 6.4.21 with the Perugia-0 tune are also shown. For FONLL and GM-VFNS the line shows the ratio 
of the central values of the theoretical cross section, while the shaded area is defined by the ratios computed from 
the upper and lower limits of the theoretical uncertainty band. 



mesons. The uncertainty on the tracking efficiency cancels completely in the ratios between production 
cross sections of mesons reconstructed from three-body decay channels (D + , D* + and D+), while a 4% 
systematic error was considered in the ratios involving the D° mesons, which are reconstructed from 
a two-particle final state. The D+/D and D+/D+ ratios were corrected for the different value of pp 
minimum-bias cross section used in and in this analysis 0. 



3 The preliminary pp minimum-bias cross section value of 62.5 mb, used in (3), was updated to 62.2 mb. 
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The predictions from FONLL (only for D°, D+ and D*+ mesons), GM-VFNS, and the PYTHIA 6.4.21 
event generator with the Perugia-0 tune are also shown 0. For all these model predictions, D mesons in 
the rapidity range \y\ < 0.5 were considered. In PYTHIA, the default configuration of the Perugia-0 tune 
for charm hadronization was used. 

The D + /D° and D* + /D° ratios are determined in PYTHIA by an input parameter, PARJ(13), that defines 
the probability that a charm or heavier meson has spin 1 . In the Perugia-0 tune, this parameter is set to 
0.54 from the measured fractions P v of heavy flavour mesons produced in vector state, see e.g. H1I29U611 . 
This setting results in an enhancement of the D + /D° and a reduction of the D* + /D° ratios with respect to 
those obtained with the default value, PARJ(13)=0.75, based on spin counting. 

The D+/D and D+/D + ratios in PYTHIA are governed by another input parameter, PARJ(2), that 
defines the s/u (s/d) quark suppression factor in the fragmentation process. In the Perugia-0 tune, PARJ(2) 
is set to 0.2, which gives rise to a reduced abundance of D+ mesons with respect to the default value of 
0.3. With this parameter adjustment, PYTHIA with the Perugia-0 tune reproduces reasonably well the 
value and pj shapes of the measured ratios involving D°, D + and D* + , while it slightly underestimates 
the abundance of D+ mesons. The fact that PYTHIA with Perugia-0 tune underestimates the strangeness 
production was already observed at the LHC in the light flavour sector 116211631 . 

In the Perugia 2011 tune 11641 . PARJ(13) is set to the same value (0.54) as in the Perugia-0 tune, while a 
lower value of the strangeness suppression factor, PARJ(2)=0. 19, is used. This results in the same values 
of the Perugia-0 tune for the D + /D° and D* + /D° ratios, and in slightly lower values for the D+ /D° and 
D+ /D+ ratios. 

The ratios of the FONLL and GM-VFNS predictions were computed assuming the perturbative uncer- 
tainty to be fully correlated among the D meson species, i.e. using the same scales for the cross sections 
at the numerator and at the denominator. Thus, the perturbative uncertainty cancels almost completely in 
the ratio, as it can be seen in Fig.|4]where, for both FONLL and GM-VFNS, the line shows the result ob- 
tained from the central values of the theoretical predictions, and the shaded area spans the region between 
the ratios computed with the upper and lower limits of the theoretical uncertainty band. The predictions 
from FONLL and GM-VFNS agree within uncertainties with the measured particle ratios. Indeed, in 
FONLL and GM-VFNS, the relative abundances of the various D meson species are not predicted by 
the theory: the fragmentation fractions /(c — > D) are taken from the experimental measurements. On 
the other hand, in both the pQCD calculations, the pj dependence of the ratios of the D meson produc- 
tion cross sections arises from the different fragmentation functions used to model the transfer of energy 
from the charm quark to a specific D meson species [28 , 65 , 66] and from the different contribution from 
decays of higher excited states. The parton fragmentation models used in the calculations provide an 
adequate description of the measured data. The measured D+/D and D+/D+ ratios do not show a signif- 
icant pj dependence within the experimental uncertainties, thus suggesting a small difference between 
the fragmentation functions of c quarks to strange and non-strange mesons. A higher statistics data sam- 
ple would be needed to conclude on a possible pj dependence of the ratios of strange to non-strange D 
meson cross sections. 

6.2 /?t -integrated D+ cross section and D meson ratios 

The visible cross section of prompt D+ mesons, obtained by integrating the pj -differential cross section 
in the measured pj range (2 < pj < 12 GeV/c), is 

a S( 2 <PT< 12GeV/c, |y| < 0.5) = 53 ± 12(stat.)^(syst.) ± 2(lumi.) ± 3(BR) jub. 

4 The ratios from the fcj-factorization model of 1161 are not shown in Fig. [4] Indeed, in this model the fragmentation 
fractions /(c -4 D) are taken from experimental measurements and the same fragmentation function is used for all D meson 
species, resulting in ratios of D meson cross sections that are independent of pj. 
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The production cross section per unit of rapidity, da /Ay, at mid-rapidity was computed by extrapolating 
the visible cross section to the full pj range. The extrapolation factor was extracted from the FONLL- 
based predictions for the D+ px-differential cross section described in Section [5] The extrapolation 
factor was taken as the ratio between the total D+ production cross section in |v| < 0.5 and the cross 
section integrated in \y\ < 0.5 and in the pj range where the experimental measurement is performed. In 
particular, the central value of the extrapolation factor was computed from the prediction based on the 
pj -differential cross section of c quarks from FONLL, the fractions f(c — > D) from ALEPH 11231 . and the 
fragmentation functions from 11571 with r = 0.1. The uncertainty on the extrapolation factor was obtained 
as a quadratic sum of the uncertainties from charm mass and perturbative scales, varied in the ranges 
described above, and from the CTEQ6.6 parton distribution functions ll67l . Furthermore, to account for 
the uncertainty on the D+ fragmentation function, the extrapolation factors and their uncertainties were 
also computed using the FONLL predictions for D°, D + and D* + mesons and the envelope of the results 
was assigned as systematic uncertainty. The resulting value for the extrapolation factor is 2.23+Qgj. The 
prompt D+ production cross section per unit of rapidity in \y\ < 0.5 is then 

dcr Di+ /dv = 118±28(stat.)+^(syst.)±4(lumi.)±6(BR)^(extr.) jUb. 

The D meson production ratios were computed from the cross sections per unit of rapidity, da /Ay. 
The corresponding values for D°, D + and D* + from were corrected to account for the updated 
value of the pp minimum-bias cross section. The systematic uncertainties on the ratios were computed 
taking into account the correlated and uncorrelated sources as described above. The resulting values 
are reported in Table @] and shown in the left-hand panel of Fig. [5] together with the results by other 
experiments that measured prompt charm production: LHCb ifTTl . e + e~ data (taken from the compilation 
in Il68l0 . and ep data in photoproduction from ZEUS 11291 and DIS from HI 11221 . The error bars are the 
quadratic sum of statistical and systematic uncertainties and do not include the uncertainty on the decay 
branching ratios, which are common to all experiments. The particle ratios for ZEUS and e + e~ were 
computed from the compilation of fragmentation fractions f(c — > D) published in [29] after updating 
the branching ratios of the considered decay channels to the most recent values 1521 . For the ZEUS 
data, the systematic uncertainties were propagated to the particle ratios by properly taking into account 
correlated and uncorrelated sources ||69l . For the HI data, the D meson ratios were computed starting 
from the unconstrained values of f(c — > D) published in ll22l . taking into account the correlated part of 
the systematic uncertainty and subtracting from the quoted 'theoretical' uncertainty the contribution due 
to the decay branching ratio IT701 . Also in this case, a correction was applied to account for the updates in 
the branching ratios of the considered decay channels. The ALICE results are compatible with the other 
measurements within uncertainties. 

The values predicted by PYTHIA 6.4.21 with the Perugia-0 tune are also shown in the figure, as well as 
those from a canonical implementation of the Statistical Hadronization Model (SHM) I7T1 . The values 
from PYTHIA were obtained by integrating the prompt D meson yields in the range \y\ < 0.5 and pj > 0. 
The SHM, which computes the hadron abundances assuming that particles originate from a hadron gas 
in thermodynamical equilibrium, provides a good description of the measured hadron yields in heavy- 
ion collisions at various energies and centralities 1721 . but it can also be applied to small systems like 
pp 117311741 and e + e~ 117511761 . The SHM results used for the present comparison were computed for 
prompt D mesons, assuming a temperature T of 164 MeV and a volume V of 30±10 fm 3 at the moment 
of hadron decoupling. The dependence on temperature of the cross section ratios considered in this 
analysis is rather small within the few MeV uncertainty on the value of T. To properly reproduce 
the yield of strange particles in small systems, such as pp and e + e~, an additional parameter, the 
fugacity 11741 . is usually introduced in the partition function to account for the deviation of strange 
particle yields from their chemical equilibrium values. For the SHM predictions reported here, a value of 
strangeness fugacity of 0.60 ± 0.04, extrapolated from the results of a fit to particle yields in pp collisions 
at yfs = 200 GeV 117711 . was used. With these parameters, the SHM provides a good description of the 
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Table 4: Ratios of the measured production cross section for prompt D mesons in pj > and \y\ < 0.5 in pp 
collisions at yfs = 1 TeV. 





Ratio ± (stat.) ± (syst.) ± (BR) 


D+/D° 
D*+/D° 
D+/D 
D+/D+ 


0.48 ± 0.07 ±0.11 ± 0.01 
0.48 ± 0.07 ± 0.08 ± 0.01 
0.23 ± 0.06 ± 0.08 ± 0.01 
0.48 ± 0.13 ± 0.17 ± 0.03 



ALICE, p T >0, |y|<0.5 

LHCb Prelim., p T >0, 2<y<4.5 

e + e\ p T >0 

ZEUS, p T >3.8 GeV/c mid-y 
H1, p T >2.5GeV/o, mid-y 
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Figure 5: Left: pj integrated ratios of D meson production cross sections compared with other experiments ffTTl 
|22][29][68]. Error bars are the quadratic sum of statistical and systematic uncertainties, without including the 
uncertainty on the BR which is common to all experiments. Right: strangeness suppression factor y s compared to 
measurements by other experiments [18, 22 . 29 . 68] . Predictions from P YTHIA 6.4.2 1 with the Perugia-0 tune and 
from a canonical implementation of the statistical hadronization model (SHM) 11711 are also shown. The gray band 
represents the uncertainty on the SHM predictions due to the uncertainty on the volume and on the strangeness 
fugacity (see text for details). 



measured ratios of D meson cross sections. 

The strangeness suppression factor for charm mesons, y s , was also evaluated. It is defined as the ratio of 
the production cross sections of charm-strange mesons (cs) to that of non-strange charm mesons (average 
of cd and cu)[§. Since all D* + and D*° mesons decay into either a D° or a D + , and all D* + decays produce 
a D+ meson Il52l . the strangeness suppression factor was computed as 



2da(D+)/dy (2) 



' da(D°)/dy + da(D+)/dy' 

The contribution to D° and D + yield from decays of excited charm-strange mesons heavier than D* + was 
neglected. 

The resulting value of y s , computed from the D+, D° and D + cross sections per unit of rapidity (da/dy), 
is 

y s = 0.31 ± 0.08(stat.) ± 0.10(syst.) ± 0.02(BR). 



5 The same symbol y s is used in the statistical hadronization model to indicate the fugacity, which, as mentioned above, is 
usually included in the partition function to account for strangeness suppression. However, the two y s are different. Indeed, in 
the statistical hadronization model, the value of the ratio between strange and non-strange charm mesons is proportional to the 
fugacity, but not equal to it, due to the different masses of the various D meson species. 
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Charm-strange meson production is suppressed by a factor « 3.3 in the fragmentation of charm quarks. 
In the right-hand panel of Fig. [51 this result is compared with the y s measurements by other experiments, 
taken from the compilation in ETTl . after updating the branching ratios of the considered decay channels 
to the values in l52l . The preliminary measurement by ATLAS [18] in pp collisions at the LHC, obtained 
using an equivalent (under the hypothesis of isospin symmetry between u and d quarks) definition of the 
strangeness suppression factor based on the cross sections of D+, D + and D* + in charm hadronization, 
is also shown. The error bars are the quadratic sum of statistical and systematic uncertainties and do 
not include the uncertainty on the decay BR. The values from PYTHIA with the Perugia-0 tune, where 
7 S corresponds to PARJ(2), and the statistical hadronization model described above are also shown for 
reference. It is also interesting to note that a similar amount of strangeness suppression was reported for 
beauty mesons by the LHCb Collaboration that measured the ratio of strange B mesons to light neutral 
B mesons, / s // d , obtaining the value 0.267+[J q2o ED- 

All the 7 S measurements, performed in different colliding systems and at different centre-of-mass ener- 
gies are compatible within experimental uncertainties. The current ALICE and ATLAS results at LHC 
energy in the central rapidity region do not allow one to conclude on a possible lifting of strangeness 
suppression with increasing collision energy. Furthermore, the D+/D (D+/D+) ratios are measured at 
the LHC both at midrapidity and at forward rapidity, thus allowing to study a possible rapidity depen- 
dence of the strangeness suppression in charm hadronization. From the comparison of the ALICE and 
LHCb results with the current experimental uncertainties (left-hand panel of Fig. [5J, it is not possible to 
draw a firm conclusion on this point. 

7 Summary 

The inclusive production cross section for prompt D+ meson has been measured in the transverse 
momentum range 2 < pj < 12 GeV/c at central rapidity in pp collisions at ^/s = 7 TeV. D+ mesons 
were reconstructed in the hadronic decay channel D+ — > <pn + with (j) — > K~K + , and charge conjugates, 
using the ALICE detector. The measured differential cross section is described within uncertainties by 
the prediction from the GM-VFNS calculation, which is based on perturbative QCD at NLO with the 
collinear factorization approach, and it is compatible with the upper side of the uncertainty band of 
calculations based on the kj -factorization approach at LO. The relative D meson production yields and 
the strangeness suppression factor, y s = 0.31 ± 0.08(stat.) ± O.lO(syst) ± 0.02(BR), agree within 
the present experimental uncertainties with those measured by other experiments for different centre- 
of-mass energies and colliding systems. More precise measurements are needed to address the possible 
energy and rapidity dependence of strangeness suppression in charm hadronization. 
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